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Abstract 


We develop a variational approximation to the entanglement entropy for 
scalar M theory in 1+1, 2+1, and 3+1 dimensions, and then examine the 
entanglement entropy as a function of the coupling. We find that in 1+1 
and 2+1 dimensions, the entanglement entropy of M theory as a function 
of coupling is monotonically decreasing and convex. While theory with 
positive bare coupling in 3+1 dimensions is thought to lead to a trivial free 
theory, we analyze a version of with infinitesimal negative bare coupling, 
an asymptotically free theory known as precarious theory, and explore 
the monotonicity and convexity of its entanglement entropy as a function of 
coupling. Within the variational approximation, the stability of precarious 
p‘^ theory is related to the sign of the first and second derivatives of the 
entanglement entropy with respect to the coupling. 
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1 Introduction 


Entanglement is a pervasive phenomenon in qnantnm physics, and has been studied 
extensively in hnite-dimensional systems using the tools of quantum information 
theory. More elusive is the entanglement structure of inhnite-dimensional systems, 
namely quantum held theories. The last twenty years has seen the development 
of a large literature on the entanglement entropy of QFT’s and has shed light on 
the entanglement between spatial degrees of freedom [D [21 a m. Currently, there 
are only a few classes of held theories for which entanglement entropy can be com¬ 
puted exactly. These include certain exactly solvable QFT’s, free held theories, and 
conformal held theories. CFT’s with holographic duals are particularly interesting 
because their entanglement entropies can be computed using the Ryu-Takayanagi 
formula [5]. 

Despite this vast body of work, the entanglement entropy of QFT’s remains 
rather mysterious. Entanglement entropies for QFT’s typically contain a combi¬ 
nation of divergent and hnite terms which depend on the parameters of the QFT. 
Often it can be difhcult to gain insight into the entanglement structure of a QFT 
by analyzing the form of its entanglement entropy. More insightful is how the en¬ 
tanglement entropy changes as a QFT undergoes dynamics such as a quench, or as 
parameters of the QFT are tuned. 

Conspicuously missing from the list of QFT’s with exactly computable entan¬ 
glement entropies are those interacting held theories which comprise the Standard 
Model. Even conventional QFT’s such as 0'^ theory and Yukawa theory have not 
had much presence in the entanglement entropy literature. While it may not be 
possible to study the exact entanglement entropy of these theories, it is possible to 
study perturbative or variational approximations. 

In this paper we embark on a program to analyze the entanglement entropies 
of more conventional interacting QFT’s, and begin with interacting scalar held the¬ 
ories and (j)'^ theory in particular. Our plan of attack is to use a variational principle 
to determine a non-perturbative approximation to the ground state of 0^ theory for 
arbitrary coupling, and then compute the entanglement entropy of the approxima¬ 
tion. Our variational ansatz will be a Gaussian wave functional, for which we can 
compute the entanglement entropy exactly. Since the ground state of a massive free 
scalar held theory is exactly a Gaussian wave functional, our approximation will be 
accurate in a neighborhood of zero coupling. We will also analyze the variational 
approximation at larger values of the coupling, for which the variational method is 
well-dehned. The validity and accuracy of the Gaussian variational approximation 
compares favorably to one-loop computations and large N approximations for theo¬ 
ries where such comparisons make sense. Our consideration of the approximation at 
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larger values of the coupling will augment our understanding of the small coupling 
case, and provide evidence for features of entanglement entropy that may hold for 
all values of the coupling. 

Once we have computed the variational approximation to the entanglement 
entropy in (j)^ theory, we will analyze how it changes as we tune the coupling. We 
will hnd that the derivatives of entanglement entropy with respect to coupling are 
often hnite and independent of regularization, and represent meaningful physical 
quantities. By studying the dependence of entanglement entropy on the coupling, 
we hnd several surprising features which suggest new insights into the nature of 
entanglement in QFT’s. 

The structure of the paper is as follows: First we discuss the variational meth¬ 
ods used to approximate the ground state and other more general states of (f)^ theory. 
Then we explain how to compute the entanglement entropy of the variational ansatz 
using the replica trick. Finally, we analyze the dependence of the entanglement en¬ 
tropy on the coupling in various dimensions for 0^ theory within the variational 
approximation. 


2 Variational Methods and the Gaussian Effective 
Potential 

In this section we introduce and review the essential aspects of the variational meth¬ 
ods we will use for the approximate calculations of entanglement entropy in inter¬ 
acting quantum held theories. The framework of variational methods that employ 
a Gaussian trial wave functional in the computation of the effective action and ef¬ 
fective potential (in particular the 2PI effective action; see Cornwall, Jackiw, and 
Tomboulis [B]) is referred to as the Gaussian Effective Potential, or GEP. The GEP 
for 0^ theory, as well as for other quantum held theories, has been studied extensively 
by Stevenson and his collaborators [3 El El do] beginning in the 1980s. That work 
built on earlier work of Barnes and Ghandour [11] and others in the 1970s, which in 
turn extended and clarihed the pioneering work of Schih [12] in the early 1960s. We 
will trace the development of the GEP and the closely related gap equation along a 
slightly diherent route, using a special case of a time-dependent variational principle 
based on work of Kerman and Koonin [13] , and applied to quantum held theory by 
Jackiw, Kerman, and others [El [181 HSl IISI HZ]- Iii considering the special case of a 
static uniform system, our results reproduce those of Stevenson et ah for scalar 0^ 
theories. 
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2.1 The Variational Wave Functional and the Variational 
Hamiltonian 


We will use the Schrodinger picture description of quantum field theory (see Jackiw 
[18] for an excellent review and additional references), where quantum held oper¬ 
ators 0(x) and their canonical momenta d(x) are time-independent, and the wave 
functional \h[0(x),f] = ((/)(x)|\[^(f)) is in general time-dependent, although we will 
not need to consider the full time-dependence in this paper. Our goal is to approxi¬ 
mately describe the quantum held theory using variational methods. We dehne the 
variational Hamiltonian functional H by 

H = {^\H\^) (2.1.1) 

where H is the Hamiltonian operator for the quantum held theory. The variational 
method consists of minimizing this variational Hamiltonian functional with respect 
to a set of variables which parameterize the wave functional T. We make the ansatz 
that the variational wave functional is a Gaussian, parameterized by a real-valued 
function (/>(x) and a real-valued symmetric kernel G(x, x') 

d'^x J d'^x' (0(x) — 0(x)) G“^(x, x') (0(x') — 0(x')) 


T[0] = iVexp 


1 
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where d is the spatial dimension. The kernel G(x, x') must be positive-dehnite 
in order for the variational wave function to be normalizable, and choosing N = 
[det(27rG)]“^/'^ assures that T has unit norm. We will frequently describe the kernel 
G(x, x') as matrix elements of an operator G acting on the space of single-particle 
wave functions. For example, in position space, G(x, x') = (x| G | x'). It is important 
to distinguish the state space of the quantum held theory where the states |T) are 
dehned and upon which the quantum held operators 0(x) and 7r(x) act, and the 
space of single-particle wave functions where the functions 0 and vr are dehned and 
upon which operators such as G and G~^ act. Using the more compact operator 
notation, we may write the variational wave functional as 


T[0] = [det(27rG)]-^/^exp 


^(0-0)G ^{ 4 >- 4 >) 


( 2 . 1 . 2 ) 


As it stands, this variational wave functional is a formal expression that will be given 
more concrete meaning through the regularization and renormalization procedures 
discussed in subsequent sections. 

Although the vacuum state of the theory corresponds to 0 = 0, other values 
of the scalar held lead to a variety of physically meaningful and interesting results. 
Considering states other than the vacuum is necessary to compute (the variational 
approximation to) the ehective action for the quantum held theory. Similar wave 
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functionals are considered in [19], where those authors were motivated by studying 
entanglement entropy in the context of free held theory in 1+1 dimensions. 

The Schrodinger picture held operators 0(x) and their canonically conjugate 
momentum operators vr(x) act on the wave functional T[0, t] in the canonical way 

0(x)T[0] = 0(x)T[0] (2.1.3) 

7r(x)T[0] = (2-1-4) 

It follows immediately that the variational parameters 0(x) are the expectation 
values of the Schrodinger picture operator in the Gaussian variational state |T) 

0(x,t) = {T(t)|0(x) |T(t)) (2.1.5) 

and we can also immediately see that G(x, x') is the two-point function or Green 
function 

G(x,x') = (T(t)l [0(x) - 0(x)] [0(x') - 0(x')] |T) (2.1.6) 


We now focus on a specihc interacting scalar quantum held theory using the 
variational method just discussed. The Hamiltonian operator for this quantum held 
theory is given by 


H 




+ \(viY + v(^) 


(2,1,7) 


where V is the scalar held potential given in terms of the bare parameters of the 
theory. We will be concerned with a quartic potential of the form 


vw 


m 


B 


+ 


A 


B 


4! 


( 2 . 1 . 8 ) 


where thb is a bare mass parameter and A^ is a bare coupling parameter. The 
coupling \b has mass dimension 2, 1, and 0 in 1+1, 2+1, and 3+1 dimensions, 
respectively. Using the Gaussian variational wave functional T and the Hamiltonian 
operator H, a straightforward calculation gives 


(T| j (i'^x-7r^(x)|\I') = ^ ~ 8 y ~ 8 I 

(2.1.9) 

For the gradient-squared term in the Hamiltonian, we hnd 

(4-1 yd-'xi(V+))"|<I>> = y'<i-*xi(V^(x))+ilVp=G (2,1,10) 


where p is the momentum operator on the space of single-particle wavefunctions, 
given by (x| p |x') = —iVx(5'^(x — x') in the position basis. Finally, the expectation 
value of the potential operator V (0) is dehned by 

F(.^,G) = (>l<(i)|W)|>l<(i)) (2.1.11) 
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where, for the quartic potential, a short computation (equivalent to applying Wick’s 
theorem) shows that 

V{^,G) =—^ + G) + — + 6 (j)'^ G + 3G ) (2.1.12) 


where we have dehned G(x) = (x| G |x) = G(x, x). Only this position-space- 
diagonal part of the Green function appears in V, because in a local quantum held 
theory all of the Schrodinger picture held operators 0(x) appearing in the Hamil¬ 
tonian operator H are evaluated at the same point in space. As in any continuum 
quantum held theory we expect G(x) to be divergent, which we will address in sub¬ 
sequent discussions of renormalization. It is worth emphasizing that the last term in 
the above expression is the square of a matrix element G = (x| G |x)^ = G(x, x)^, 
not the matrix element of the square of the operator G, which we would write as 
(x| G^ |x) = G^(x, x). Combining the results above, the full variational Hamiltonian 


IS 


H = Tr 






-I- / d‘^x 


-(V^f + V(^,G) 


(2.1.13) 


The equations for the variational parameters, for an arbitrary potential function V, 
follow directly from minimizing this functional with respect to 0 and G 


dcj) 


(2.1.14) 

5H 

4 

(2.1.15) 


where the last equation is written in the condensed operator form. The operator u 
is dehned by its kernel in position space as 

cj^(x, x) = (—Vx + d^(x — x') (2.1.16) 


where the quantity G) is dehned, for an arbitrary potential function V, by 


m^ = 2^ (2.1.17) 

dG ^ ^ 

and for the quartic potential in particular, it is 

= ‘^^ = ^B + ^iP + G) (2.1.18) 

oG 2 

The quantity m will play a signihcant role in subsequent developments. In addition 
to depending on the bare mass and bare coupling A^, m is a function of 0(x) 
and G(x), so it may depend on position in general. In the case of a uniform system, 
m will be a constant parameter. We now turn our attention to such static uniform 
systems. 
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2.2 The Gap Equation and the Gaussian Effective Potential 


If we now restrict ourselves to spatially uniform systems, which are translation 
invariant, then 0 does not depend on position and = 0. Additionally, the Green 
function G(x, x') will only depend on the difference x — x'. Because of translation 
invariance, the analysis simplifies considerably in momentum space. Defining the 
Fourier transforms 

and using the definition of of cu, its Fourier transform takes the simple form 

u{p,m) = (2.2.3) 


( 2 . 2 . 1 ) 

( 2 . 2 . 2 ) 


where we now indicate explicitly the dependence m. The second of the two non¬ 
trivial Hamilton equations, for G~^ in terms of is now 

^ G(p, = a;(p, m)^ = -I-(2.2.4) 

This important equation, together with 

rri^ = rn% + + G{m)] (2.2.5) 


is called the gap equation. The expression for G, now given in terms of m, is 


G(m) = G(x, x) = 


d'^p 

( 2 ^ 


G(p,m) = 


d'^p 


(27r)'^2a)(p, m) 


( 2 . 2 . 6 ) 


The gap equation gives the functional form of the Green function G in the static 
uniform case and an equation connecting the mean field <p to the mass parameter 
m. We will frequently refer to just this relationship between 0 and m in (I2.2.5j) 
as the gap equation, although the full meaning is as just discussed. The Green 
function describes a particle in the interacting quantum field theory propagating 
with a mass m. We have arrived at this interpretation of the Green function G and 
the mass m via the variational principle and the Gaussian form of the variational 
wave functional; it is not a separate assumption. 

The variational Hamiltonian simplihes in the uniform case to 




+ 


id = Tr 


d'^x!/((/., G) 


(2.2.7) 











Computing the trace in the momentum space representation, the variational Hamil¬ 
tonian density "H, dehned hy ff = f is given by 




{27iY 


iG(p)-‘ + lp^G(p) 


+ V(lG) 


( 2 . 2 . 8 ) 


If we look back through the derivation of this expression, we see that the first term, 
with G~^, comes from the expectation value of the term of the Hamiltonian 
operator; the second term, with p^G, comes from the expectation value of the (V0)^ 
term of the Hamiltonian operator; and the third term V comes from the expectation 
value of the potential V (0) term of the Hamiltonian operator. The sum of the first 
two terms is less divergent than either one separately, which is a consequence of the 
Lorentz invariance of the theory. We now use the relationship of the momentum 
space Green function G(p,m) to a;(p, m) via the gap equation, so the momentum 
space integrals may be written as 


d'^p 


^G'(p,m) ^ + ^p2G(p,m) 


d'^p 


^G(p,m) ^ - ^G(p,m) 


d'^p 


2 m- 


m 


J (27r)'^ 2ci;(p, m) 

where, following [9], we have dehned the integrals 


= Ii{m) - —/o(m) 


lN(m) = 


d^P X2JV 1 

-—-1 -;--a;(p,m) =- 

(27r)‘^ 2a;(p, m) 2 


/p2 + ^2)7V-i ^2.2.9) 


(27r 


We note in particular that G{m) = Io{m) is the integral appearing in the gap 
equation and the V part of the variational Hamiltonian. Most of these integrals will 
be divergent in the cases we will be studying, and the renormalization prescriptions 
we are about to discuss will address this in detail. 

The variational Hamiltonian density "H, with m given as a function of the mean 
held 0 via the gap equation, is known as the Gaussian Effective Potential (GEP), 
and is given by 

2 

_ Tf) _ 

VGi?(0) = h{m) - —/o(m) V (2.2.10) 

and V for the quartic potential was given earlier by 

^ ^ (0"^ 60^/o(m) + 3 Io{mY) (2.2.11) 


where we have used G{m) = Io{m). 

We now turn to the renormalization of the theory and a physical interpretation 
of the formal expressions for the gap equation and the GEP. 
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2.3 Renormalization of the Gaussian Effective Potential 


The vacuum state, with (/> = 0, corresponds to a minimum of Vce{<P) 

dVcE 


d(p 


= 0 


(2.3.1) 


b=0 


Following Stevenson [Hj as well as Barnes and Ghandour P, we define the renor¬ 
malized mass ruR 

2 d'^VcE 

mR = 


dcj)'^ 


(2.3.2) 


A=0 


and the renormalized coupling Xr 


— 


d^V( 


GE 


d(j)^ 


(2.3.3) 


6=0 


Computing these derivatives is routine, and the results are quite simple when the 
derivatives are evaluated in the vacuum state 0 = 0. For the renormalized mass niR 
we have 


A 


'^R = '^B+ '-Y ^o("io) = ml 


(2.3.4) 


The second equality above follows because mo satishes the gap equation m^ = 
m^ + ^ do{mo) with 0 = 0. So, the renormalized mass niR is identified with the 
vacuum solution of the gap equation: rriR = mg. The renormalized mass sets the 
scale for all dimensionful quantities in the renormalized theory, and it will often be 
convenient to work in units such that niR = 1. For the renormalized coupling Xr, 
we have 

X _x 1 - AB/-i(mij)/2 

■^R — :r—w— j —7-^ 

1 -h XBl-i{mR)/A 

With these renormalization prescriptions, we can now derive the renormalized forms 
of the gap equation and the GEP in various dimensions. 

The detailed discussion of the 1-1-1 and 2-1-1 dimensional cases are quite similar. 
In these dimensions there are ultraviolet divergences that require mass renormaliza¬ 
tion, but coupling renormalization is merely a hnite coupling reparameterization, 
and in these lower dimensions we will express results in terms of the original finite 
bare coupling Xr- We refer the reader to the papers of Stevenson and collaborators 
for more detail, but we will review and summarize the most important results here. 
In 1-1-1 and 2-|-l dimensions, we will make use of the following identities for the 
divergent integrals Ii{m) and Io{m): 


/i(m) - h{mR) = -{w? - m\) hirnR) + Fi(m) 

/o(m) - Io{mR) = Fo(m) 


(2.3.6) 

(2.3.7) 
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where the finite quantities FQ{m) and Fi{m) are 


Fi(m) 

1 / 2i 2 

= - —(m log - m + I) 

Svr 

I + I dimensions 

(2.3.8) 

Fi(m) 

I 

{m — mB)‘^{2m + mji) 

2 + I dimensions 

(2.3.9) 

247r 


Fo{m) 

I 

I + I dimensions 

(2.3.10) 


Fo{m) 

I 

I 

2 + I dimensions 

(2.3.ll) 


The renormalization of the gap equation proceeds in essentially the same way in 1+1 
and 2+1 dimensions. Starting with the bare gap equation and using the definition 
of the renormalized mass, we have 

= m| + ^/o(m) + (2.3.12) 

= "^1 + + ^(/o(m) - loimn)) + ^0^ (2.3.13) 

= "i/? + ^^o(m) + (2.3.14) 


As a check, note that when 0 takes its vacuum value ^ = 0, the mass parameter m 
also takes its vacuum value m = mo = m^, since FoirriR) = 0. So, the renormalized 
gap equations for each dimension are 


72 1 , 2,2 2 \ 

6 = — log - 2 - + T—(m^ - m^) 
dvr \b 

1 2 

= — (m — m^) + -— {m? — mj^) 
dvr Xb 


1 + 1 dimensions (2.3.15) 

2 + 1 dimensions (2.3.16) 


In either case, the renormalized gap equation is of the form 

2 

As 


= —Fnim) + -—(m^ — 


R) 


(2.3.17) 


We see the similarity of the form of the gap equations across different dimensions. 
Only the first term has a different form in different dimensions and has no de¬ 
pendence on the form of the potential in the original Hamiltonian. The last term 
contains the coupling \b and depends on the detailed form of the potential. As we 
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will see, a very similar result holds in 3+1 dimensions, where the gap equation will 
be expressed in terms of the renormalized coupling instead of the bare coupling. 

The renormalization of the Gaussian Effective Potential Vge also proceeds in 
essentially the same way in 1+1 and 2+1 dimensions. Consider the hrst two terms 
oiVcE-. 


2 2 
TTl Tfl 

h{rn)—= h{mR) + /i(m) - h{mR) - —Io{m) 

1 TTh^ 

= Ii{mR) + -(m^ - m\)lQ{mR) + Fi(m)- 

2 2 

= hiniR) - ^ h{mR) + Fi{m) - ^(Jo(m) - /o(m^)) 

= Ii{mR) - ^ h{niR) + Fi(m) - ^Fo(m) 

The hrst two terms contribute an additive divergent constant to Vgb, but we will 
require the energy density to be zero in the vacuum state, using the freedom to add a 
constant to the original energy density to make this so. Using the above expressions 
for Fq and Fi, the hnite parts are quite simple: 


, , W? ^ , 1/9 9 X 

= — (m - m^) 

Tfl^ 1 


1 + 1 dimensions (2.3.18) 

2 + 1 dimensions (2.3.19) 


Turning to the third and hnal term of Vge, namely V, a useful form is given by 
(recall that G(m) = Io(m)) 


-2 A 


a) = ^(f-‘ + G(m)) + + G(m) + 3 a(mf) 


74 , 1- / 2 , \\ 72 , 1- / 2 , ^R 


m 


R 


4! 

A_B 74 1 


+ ^i^R + ^ ^ G{m)) 


4! 


H—(m^ 
2^ 


Ab 72\ 72 I 1/2 A_Bt2\2 

+^(™ -^4>) 


2A 


B 


_ 74 , 

12 ^ ^ 2 Xb 2 \b 

where we have once again used the bare gap equation. This expression fl2.3.20p 
holds in any dimension. Note that the last term above is also an additive divergent 
constant. Combining the results above for the two parts of Vgb, and choosing the 
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overall additive constant such that Vge = 0 in the vacuum state characterized by 
0 = 0 and m = mo = m^, we have the hnal result for the renormalized GEP in each 
dimension 


Vge{(I>) = - ml) + 


m 


R) 




m 


R) 


2A 


[m — m 


B 


R) 


12 ^ 

12 ^ 


1 + 1 dimensions 

(2.3.20) 

2 + 1 dimensions 

(2.3.21) 


As already noted for the gap equations, this way of writing things clearly shows the 
similarity of the GEP across different dimensions. Only the hrst term has a different 
form in different dimensions; it arises purely from expectation value of the operators 
TT^ and (V0)^ terms in the Hamiltonian operator H, and has no dependence on the 
form of the potential operator V (0). The last two terms, which contain the coupling 
Xb, depend only on the choice of potential. Of course, the dependence of m on 0 
is given by the appropriate gap equation in each dimension, and these have similar 
but not identical forms. 

In 3+1 dimensions, we will restrict our attention to a version of 0^ theory 
which is non-trivial in 3+1 dimensions - meaning it is physically well-dehned and 
is not a free held theory. This is the theory which Stevenson [9] has given the name 
precarious 0^^ theory. In this theory, the bare coupling is taken to be negative and 
arbitrarily close to zero in a limit in which the ultraviolet cutoff goes to inhnity. 
This theory, while unstable, is in a very specihc sense arbitrarily metastable as the 
ultraviolet cutoff goes to inhnity, and constitutes a physically sensible ehective held 
theory. We do not attempt to review precarious 0^ theory and refer the reader to the 
very clear discussion in Stevenson [9]. As noted there, and observed much earlier by 
Symanzik [20] and Parisi [2T], this version of 0^ theory has the intriguing property 
of being asymptotically free. We do not consider the “standard” version of 0^ theory 
with positive bare coupling because it is trivial using the variational methods of the 
GEP - meaning it is a free held theory with the renormalized coupling driven to zero 
from above as the ultraviolet cutoh goes to inhnity. This is in harmony with other 
“triviality” results for standard 0^ theory in 3+1 (and higher) spacetime dimensions. 
According to Stevenson and Tarrach [221 123] , in addition to the precarious version 
there is one other physically sensible version of 0^ theory in 3+1 dimensions which 
has a positive bare coupling but dihers from the standard theory. They call this 
theory autonomous 0^ theory. We will not consider this version in this paper. 

In 3+1 dimensions the ultraviolet divergences require coupling renormalization 
as well as mass renormalization. The prescription for the renormalized coupling in 
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any dimension, as formulated earlier, is 


d4>‘^ ^0=0 ^ ^ 


(2.3.22) 


In the precarious (p'^ theory, since the bare coupling As is negative and driven to 
zero as the cutoff goes to inhnity, the coupling renormalization is given by 


XbI _i(ms) 
4 


-1-3 


4 

Xr /_i(ms) 


+ ■■■ 


(2.3.23) 


where we have dropped terms of order J_i(ms)“^ or smaller. The following iden¬ 
tities for the integrals /Ar(m) in 3-1-1 dimensions are useful for the renormalization 
procedure: 

Ji(m)-/i(ms) = ^ (m^ - m%) /o(ms) - ^(m^ - /_i(ms) + Fi(m) (2.3.24) 

Z o 

Io{m) - /o(ms) = I-i{rnR) + Fo(m) (2.3.25) 

I-iim) - /„i(ms) = F_i(m) (2.3.26) 


where the hnite parts are 


1 ?77^ 

-3("i^-m|)^-2m|(m2-m^)) (2.3.27) 

1287r^ m\ 


. \ 1.9 2 2i \ 

Fo m) =-—^(m -m^-m log—) 

^ . 1 


(2.3.28) 

(2.3.29) 


Using these identities and the coupling renormalization, it is now straightforward to 
derive the renormalized gap equation, which is 


1.0 2 2 1 ^ / 2 2 \ 
(m -m^-m log—) + —(m - m^) 
lOTT^ Xr 


(2.3.30) 


Following similar derivations already discussed for lower dimensions, a straightfor¬ 
ward calculation shows that the renormalized GEP is 


m 


1287r2 


m 


R 


Vge{,^) = 77 ^ log “ ^r) 


(2.3.31) 


It is worth noting that there is a non-trivial cancellation of logarithmic divergences 
between the kinetic and potential pieces of the Hamiltonian, associated with the 
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coupling renormalization. This is the key to the renormalizability of the theory in 
3+1 dimensions. It does not occur and is not required in lower dimensions. 

We now express all dimensionful quantities in units of the renormalized mass 
rriR, so that = 1. It is convenient to dehne the normalized coupling a = 
and the negative normalized inverse coupling k = —1/a, and a re-scaled version of 
the scalar held by </? = 47r(/ so that (p'^ = 167r^(/^. Then the renormalized gap 
equation takes the form 

p? = ((k — 1)(1 — m?) — logm^) (2.3.32) 

and the renormalized GEP is 

Vge = - 1)(1 - - m^logm^) (2.3.33) 

Following [9], for each value of the coupling we dehne the critical value pc as the 
largest value of the scalar held where the gap equation has a solution. These are 
given by 


pI = K — \ + exp(— k) (2.3.34) 

= exp(— k) (2.3.35) 


The break value pb is dehned as the largest value of 0 where the solution of the 
gap equation gives the global minimum of the GEP. When p > pb, the GEP is a 
constant Vge{p) = VGE^Pb), corresponding to the minimizing value of m = 0. 


pI = K - 1 + - exp 



ml = exp 



(2.3.36) 

(2.3.37) 


It is always the case that pb < Pc and mb > me- The region between these val¬ 
ues corresponds to a metastable region of the GEP, where the solution to the gap 
equation corresponds to a local but not global minimum. As discussed in detail in 
[9], when the coupling is strong enough, specihcally when k <1/2, the break value 
Pb goes to zero and the GEP becomes flat for all values of p, corresponding to an 
interacting theory with a massless particle. Keeping track of the critical values and 
break values will be important in the discussion of our results for the entanglement 
entropy in 3+1 dimensions. 

We need to make some hnal observations before we move on to entanglement 
entropy, concerning the sign of the coupling in each dimension. Once again, Steven¬ 
son [9] provides a clear discussion of these issues, and we merely summarize his 
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findings. In 1+1 and 2+1 dimensions, the concern is over infrared divergences. In 
these cases it is essential that the bare coupling be non-negative: > 0. Other¬ 

wise, the Gaussian Effective Potential is either completely ill-dehned (in 1+1) or is 
unbounded below (in 2+1). In 3+1 dimensions for precarious (j)'^ theory the concern 
is not infrared divergences and a GEP unbounded below, but ultraviolet divergences 
that require coupling renormalization. However, something even more interesting 
takes place: there is a duality relation that connects positive renormalized couplings 
to negative renormalized couplings, along with a rescaling of the mass parameter m 
and the scalar held 0, and with an additive constant shift of the GEP. The explicit 
formulas dehning the duality transformation are 

ml=m‘^/F = fi;* = K + logF (2.3.38) 

where the asterisk denotes the dual values, and the scale factor F is the solution of 

(F-1)(1-k) = FlogF (2.3.39) 

The entire range of positive couplings a > 0 (or equivalently negative inverse cou¬ 
plings K < 0) gets mapped to a range of negative couplings a < — 1 (or equivalently a 
range of positive inverse couplings 0 < k < 1) under duality, so the duality mapping 
is one-to-one but not onto. Note the somewhat counter-intuitive results that k = 0 
is the self-dual point, indicating that the two formulations with inhnitely negative 
coupling and inhnitely positive coupling are dual to each other, and that arbitrar¬ 
ily weak positive coupling is not dual to arbitrarily weak negative coupling. This 
duality relation means we need only consider negative values of the renormalized 
coupling a with no loss of generality. In later sections on entanglement entropy we 
will express our results in terms of the negative renormalized couplings, although 
we will briehy discuss how some of these results look for positive couplings which 
are dual to the negative couplings for purposes of comparison with results in lower 
dimensions. 

On the following pages we plot the solutions of the renormalized gap equation 
and the renormalized Gaussian Ehective Potential in 1+1, 2+1, and 3+1 dimensions, 
for diherent values of the coupling. We also plot the duality relations for coupling 
and scale. We emphasize that all dimensionful quantities in these plots are measured 
in units of the renormalized mass scale mn, which sets the scale for the entire theory. 
In these units, mu = 1. 
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Gap Equation Solutions in 1+1 Dimensions 
minimum lambdaB/4! = 1 in red; maximum lambdaB/4! = 5 in black 



Figure 1: Solutions of the renormalized gap equation in 1+1 Dimensions, with a 
range of values for the bare coupling Xb- The bottom curve (red in color plot) 
corresponds to the weakest coupling; the top curve (black in color plot) corresponds 
to the strongest coupling. 
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Mass Parameter m 


Gap Equation Solutions in 2+1 Dimensions 
minimum lambdaB/4! = 1 in red; maximum lambdaB/4! = 5 in black 



Figure 2: Solutions of the renormalized gap equation in 2+1 Dimensions, with a 
range of values for the bare coupling A^. The bottom curve (red in color plot) 
corresponds to the weakest coupling; the top curve (black in color plot) corresponds 
to the strongest coupling. 
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Gap Equation Solutions in 3+1 Dimensions 
minimum kappa = 0.8 in red; maximum kappa = 1.2 in black 



Figure 3: Solutions of the renormalized gap equation in 3+1 Dimensions, with a 
range of values for the inverse coupling k. The innermost curve (red in color plot) 
corresponds to the strongest coupling (most negative a or equivalently smallest k); 
The outermost curve (black in color plot) corresponds to the weakest coupling (least 
negative a or equivalently largest k). 
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Gaussian Effective Potential 


Gaussian Effective Potentiai in 1+1 Dimensions 



Scalar Field 


Figure 4: Renormalized Gaussian Effective Potential in 1+1 Dimensions, with a 
range of values for the bare coupling A^. The top curve (red in color plot) corre¬ 
sponds to the weakest coupling; the bottom curve (black in color plot) corresponds 
to the strongest coupling. 
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Gaussian Effective Potentiai in 2+1 Dimensions 
minimum lambdaB/4! = 1 in red; maximum iambdaB/4! = 5 in biack 



Figure 5: Renormalized Gaussian Effective Potential in 2+1 Dimensions, with a 
range of values for the bare coupling Xb- The top curve (red in color plot) corre¬ 
sponds to the weakest coupling; the bottom curve (black in color plot) corresponds 
to the strongest coupling. 
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Gaussian Effective Potentiai and Scaiar Fieid in 3+1 Dimensions 
minimum kappa = 0.8 in red; maximum kappa = 1.2 in biack 



Figure 6: Renormalized Gaussian Effective Potential in 3+1 Dimensions, with a 
range of inverse negative couplings k. The bottom curve (red in color plot) corre¬ 
sponds to the strongest coupling (most negative a or equivalently smallest k); The 
top curve (black in color plot) corresponds to the weakest coupling (least negative 
a or equivalently largest k). 
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coupling k 



Figure 7: Duality relation for precarious negative inverse coupling k > 0 as a 
function of dual inverse coupling k* < 0. The duality transformations are given in 
equations (12.3.381) and (12.3.391) . 
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duality scale factor 



Figure 8: Duality relation for precarious 0^: duality scale factor F as a function of 
dual inverse coupling k* < 0. The duality transformations are given in equations 


fl2.:-i.:-i8D and (12X39]). 
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3 Entanglement Entropy from the Replica Trick 


Since our variational ansatz for theory takes the same form as the ground state of 
a massive free scalar held, we would like to know how to compute the entanglement 
entropy of a massive free scalar held. We will use the replica trick, which pro¬ 
vides an efficient way of computing the von Neumann entropy of a density matrix. 
Let p denote an un-normalized density matrix, and let p denote the corresponding 
normalized density matrix p := p/Tr(p). Then the replica trick [26] gives us 

lim (—^ -I- 1 ) logfTrp") = — Trfplogp) (3.0.1) 

\ an J 


equation fId.O.ip is particularly useful since it normalizes our density matrix of in¬ 
terest, p. This allows us not to worry about the normalization of our density matrix 
prior to applying the replica trick. 


With the help of the replica trick, computing the entanglement entropy is con¬ 
ceptually straightforward but technically difficult. Given a wave functional T[0(x)] 
of our quantum held theory, we start by tracing out the degrees of freedom in a 
spacetime region V from the density matrix T[0(x)]T*[0(x)] to obtain the reduced 
density matrix p, which need not be normalized. Then we compute Tr p"", and use 
the replica trick to obtain the desired entanglement entropy. 


In our case, we consider the Gaussian wave functional discussed in detail earlier 


= [det(27rG)] 


-1/4 


exp 




Although we are actually interested in wave functionals which take the form 
T[(/)(x)] = \h[0(x)-|-0], letting p be the reduced density matrix of T[(/)(x)]T*[0(x)] = 
T[(/)(x) — 0]T*[0(x) — 0], we see that Tr p" = Tr p^ by a change of variables. There¬ 
fore, it suffices to consider the entanglement entropy corresponding to In the 

context of free held theories, the entanglement entropy for such Gaussian states 
has been considered by a number of authors, originating in the work of Sorkin and 
collaborators pTl [25] . 


For this paper, we will restrict ourselves to the case in which V is the half¬ 
space. When V is the half-space, it can be shown that Tr p” is proportional to the 
partition function of a massive free scalar held living on an n-sheeted covering 
of the Euclidean plane, which is a cone of dehcit angle 5 = 27r(l — n). This dehcit 
angle is non-positive because each power of the density matrix adds another sheet 
which contributes 27r radians of angle “surplus.” Since ^ = —2%^, when V is the 
half-space equation fl3.0.ip can be written as [26l |28] 


lim 

5^.0 


d 


^Og Zs 


Tr(plog p) 


(3.0.2) 
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Using heat kernel methods, it is possible to compute logZ^, which will give us 
the leading area law contributions to the entanglement entropy [271 [281 [29]. Using 
equation fl3.0.2p we hnd that 

S = — ,,, d + 1 dimensions (3.0.3) 

12(47r)U-i)/2 sU+i)/2 ^ ^ 


where e is the length scale which serves as our UV cutoff, S is the boundary of V, 
and A is the area of S. The expansion of equation fl3.0.3p contains divergent terms 
due to the UV cutoff. In this paper, we will be interested in the cases d = 1,2, 3. 
For d = 1, we have 

S=-hE.(-mV) 

= (log(™^e^) + 7 b) + 0(m?iL‘) 

= —— (loge^ + logm^) + 0{w?e^) 1 + 1 dimensions (3.0.4) 

where Ei( 2 ;) is the exponential integral and is the Euler-Mascheroni constant, 
and without loss of generality we have rescaled e to e = . In d = 1 dimensions 

only, we have that A = 1. For d = 2, 


s = -l- 

12 I vTf 


m + merf(me) 


A 


— m] + 0{rn?e) 


12 V0re 

where eTf{z) is the error function. Finally, for d = 3, 


2 + 1 dimensions 


(3.0.5) 


F = 


A 


+ Ei (—m^e^) 


A 


487r V 


A 


— + log log (y^; — 1) ) + 0{me) 


487r V 


37 b 


+ w? log + m? log m? — wA ) + 0{me) 


1 dimensions 
(3.0.6) 


where, as in 1+1 dimensions, we have again rescaled e to £ = 

We note that there are other sub-leading area terms which can be obtained 
using more sophisticated methods [30], but they will not affect our analysis. In 


26 














particular, when we consider derivatives of the entanglement entropy with respect 
to the coupling, these other sub-leading area terms will either vanish in the e —?■ 0 
limit, or if they diverge will be sub-leading to other divergences. 

One may notice that the divergences in S' in d spatial dimensions take the same 
form as the divergences of the effective potential for the massive free scalar held in 
d spatial dimensions. Indeed, this is no coincidence. Noting that Wg = — log Zg is 
the effective action for a free massive scalar held with dehcit angle 6, we can rewrite 
equation fl3.0.2p as 


lim (-27r^ - 1 ] = - Tr(plogp) (3.0.7) 

<5->-o \ do / 

Since the ehective action contains the ehective potential and thus the divergences of 
the ehective potential, it is unsurprising that the entanglement entropy also contains 
divergences of the same form. In fact, the entanglement entropy takes the following 
simple unifying form in all dimensions: 

S/A = I G{m) = ^/o(m) (3.0.8) 

where we recall that G{m) = G(x, x) is the divergent diagonal piece of the Green 
function in position space. Let us dehne AS by the hnite part of the entanglement 
entropy S plus an additive constant which makes AS equal to zero in the free held 
vacuum. Then in all dimensions we have 

AS/A=^Fo{m) (3.0.9) 

For our variational approximation, we will simply replace m in the above equations 
with the variational mass. 


4 Calculations of Entanglement Entropy 

4.1 Entanglement Entropy Per Unit Area 

In this section, we will compute the entanglement entropy of 0“^ theory within the 
variational approximation, and examine how the entanglement entropy depends on 
the coupling. In a preceding section on the variational method and Gaussian Ehec¬ 
tive Potential, we found that the renormalized gap equations in spatial dimensions 
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(i = 1, 2, 3 take very similar forms, which we repeat here: 


— log —^ + -—(m — mr/j 1 + 1 dimensions 

471 Xb 

(4.1.1) 

1 2 

— (m — mji) + -—(m^ — m^) 2 + 1 dimensions 

47r Xb 

(4.1.2) 

- - (ttib — log -^) — — rn?) 3 + 1 dimensions 

IGvr^ Xr 

(4.1.3) 

Also, in the previons section on the formalism of entanglement entropy, we fonnd 
that expressions for AS in spatial dimensions d = 1, 2, 3 also take very similar forms, 
expressed in all cases by AS/A = ^Fo{m): 

A5M = -llog^ 

12 mj^ 

AS/A = (m - rriR) 

^ log ^ ) 


To pnt all of these eqnations in the exactly the same form is now jnst a qnestion 
of normalization, so we dehne the normalized dimensionless conpling a and the 
normalized dimensionless scalar held for each dimension as 


a = 

As 

ip^ = 47r0^ 

1 + 1 dimensions 

(4.1.7) 

87rm|, 

a = 

As 

= 47r0^ 

2 + 1 dimensions 

(4.1.8) 

STirriR 



a = 

As 

= 167r^0^ 

3 + 1 dimensions 

(4.1.9) 

and the normalized 

dimensionless qnantities appearing in the gap eqnation 



/ = 

log rn? 

1 + 1 dimensions 

(4.1.10) 


/ = 

m — 1 

2 + 1 dimensions 

(4.1.11) 


/ = 

— (1 — logm^) 

3 + 1 dimensions 

(4.1.12) 


1 + 1 dimensions (4-1.4) 


2 + 1 dimensions (4.1.5) 


3 + 1 dimensions (4.1.6) 








where m and other dimensionful parameters are now expressed in units of the renor¬ 
malized mass scale rriR, so in effect we have chosen units such that rriR = 1. Finally, 
define the normalized entanglement entropy per unit area a by 


a= 12AF/H 

1 + 

a = 12 AS/A 

2 + 

a = 487r AS/A 

3 + 


1 dimensions (4.1.13) 

1 dimensions (4.1.14) 

1 dimensions (4.1.15) 


The term “area” is used to refer to the measure of the surface of spatial co-dimension 
1 in all cases, and as with other dimensionful quantities it is measured in units of 
the appropriate inverse power of the renormalized mass scale In terms of these 
quantities, the dimensionless gap equation and the normalized entanglement entropy 


per unit area take the same simple form in any dimension: 

= firri^) H- — 1) (4.1.16) 

a 

a = -/(m^) (4.1.17) 

We note the ranges of the relevant quantities in various dimensions are 

a > 0 > 1 >0 1-1-1 dimensions (4.1.18) 

a > 0 > 1 /(’^^) > 0 2+ 1 dimensions (4.1.19) 

a < 0 0 < m? <1 — 1 < fim^) < 0 3-|- 1 dimensions (4.1.20) 


The system of equations given by (I4.1.16p and (I4.1.17p allows us to determine the 
variational approximation to the entanglement entropy for fixed values of the cou¬ 
pling and vacuum expectation value p. In particular, for a hxed value of the coupling 
and a fixed value of p, we solve the gap equation in Eqn. fl4.1.16p for the varia¬ 
tional mass m, and then plug m into equation fl4.1.17p to obtain the normalized 
entanglement entropy per unit area a. 

Before delving into a full discussion of the dependence of entanglement en¬ 
tropy on the coupling, we briefly point out a connection between our variational 
calculation of entanglement entropy in 0“^ theory and a perturbative calculation. 
For d = 3, 0'^ theory is widely believed to be a sick theory for finite positive cou¬ 
pling in a non-perturbative sense. However, one can make sense of the theory for 
inhnitesimal positive coupling, and treat the theory perturbatively. To compute the 
entanglement entropy in the vacuum of such a theory, we consider equation (I4.1.16p 
for d = 3 and </? = 0 in the a —)■ 0 limit, and find that = 1. Temporarily restor¬ 
ing the renormalized mass scale we find that the solution to the gap equation 
in the a —)■ 0 limit is m = rriR. Then the approximate normalized entanglement 
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entropy per unit area is a = = 0, which agrees exactly with a perturba¬ 

tive calculation of Hertzberg to hrst order in perturbation theory [31]. This result 
demonstrates that our variational framework can reproduce perturbative results in 
the small coupling limit. However, we note that in the case of (j)"^ theory with in- 
hnitesimal positive coupling, we learn nothing about the explicit dependence of the 
entanglement entropy on the coupling since a = = 0 does not depend on 

the normalized coupling a. 


4.2 Dependence of Entanglement Entropy on the Coupling 


Now we will analyze the dependence of the normalized entanglement entropy per unit 
area a on the coupling. The results presented in this section demonstrate the utility 
of the variational techniques, and explore the relationship between entanglement 
and coupling. 


In order to determine a as a function of coupling, there are two approaches. 
The hrst approach is to solve the gap equation (14.1.161) at hxed coupling and hxed 
If for the variational mass m, and then plug m into equation (14.1.171) to obtain a. 
In fact, the gap equation (14.1.161) can be solved exactly in terms of Lambert W 
functions for d = 1 and d = 3, and by the quadratic formula for d = 2. However, 
the form of the corresponding normalized entanglement entropy per unit area a is 
not particularly enlightening since it is difficult to parse, although we can graph it 
numerically and examine the plot. The normalized entanglement entropy per unit 
area a is plotted in each dimension on the following pages. 


While we can guess some trends from the plots, it is best to be more quantita¬ 
tive. A second approach to characterizing the dependence of a on the coupling is to 
compute (1^)^ and implicitly using the gap equation, and then examine the 

behavior of these derivatives. Following this approach, let us hrst examine 
First we note that 




dm?' 


\ 


da 


(4.2.1) 


While we can easily compute = df /dm? using the explicit form of / in 


each dimension, we need to use the gap equation to compute ■ Considering 

m as a function of the scalar held and coupling given by the gap equation, and 
diherentiating the dimensionless gap equation with ip held hxed leads to the general 
result 



— 1 / 1 
a^ \a 



(4.2.2) 
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Specifically, in each dimension, 


dm'^ 

— 1 

1 + 1 dimensions 

(4.2.3) 

da 

a (1 + a/m^) 

dm^ 

— 1 

2 + 1 dimensions 

(4.2.4) 

da 

«(1 + a/2m) 

dm? 

m? — 1 

3 + 1 dimensions 

(4.2.5) 

da 

a (1 + cc logm^) 


Calculation of the hrst derivative of the normalized entanglement entropy per unit 
area a is now simple: 




da j, 


(m^ — 1) 


a^ 


1 + 


af'{m? 


-1 


(4.2.6) 


Using the ranges of the relevant quantities noted above, we find 


da 




< 0 


da 


‘p 


< 0 


da 




< 0 


da 




> 0 


1 + 1 dimensions 

(4.2.7) 

2 + 1 dimensions 

(4.2.8) 

3 + 1 dimensions, w? > for a < 0 

(4.2.9) 

3 + 1 dimensions, m? < for a < 0 

(4.2.10) 


We also note the results for positive coupling a in 3+1 dimensions, which are related 
to the negative coupling results by the duality relation we discussed previously 


® *• 


3 + 1 dimensions, for a > 0 




( 7 — I <0 3 + 1 dimensions, for « > 0 


(4.2.11) 

(4.2.12) 


We will now examine the results in equations (I4.2.7p through (I4.2.10p . First, 
for equations fl4.2.7p and fl4.2.8p . we note that 


_ U /dAS\ _ U 
da)^ A\da)^ A [daj 


1+1 and 2+1 dimensions (4.2.13) 
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In the above equation, we can replace AS* with S since the divergent pieces of S 
have no dependence on the coupling for d = 1 and d = 2. Thus, while S is divergent 
and regularization dependent, is finite and independent of regularization, and 

so represents a meaningful physical quantity. On the other hand, for precarious 0'^ 
theory where d = 3, is not hnite and so is not proportional to However, 

the divergent term in (|f)^ is proportional to log(e) and so in light of equation 
fl4.2.5p . for small e we hnd that (ff)^ has the same sign as (|^)^- So to analyze the 
monotonicity of S for precarious 0^ theory, it suffices to analyze the monotonicity 
of a. 

In equations fl4.2.7p and fl4.2.8l) . we see that for d = 1 and d = 2 the variational 
approximation to the entanglement entropy is monotonically decreasing with respect 
to coupling for fixed (p. This result is surprising, since one might have expected the 
entanglement entropy to increase as the coupling is increased. We can trust our 
result in a neighborhood of zero coupling for which the variational approximation 
is accurate. Thus, a conservative interpretation of equations fl4.2.7p and fl4.2.8p is 
that for d = 1 and d = 2 the entanglement entropy is monotonically decreasing in a 
neighborhood of zero coupling. More ambitiously, we might take the monotonicity of 
the variational approximations of the entanglement entropies as an indication that 
the true entanglement entropies for d = 1 and d = 2 are monotonically decreasing 
for all values of the coupling. This amounts to little more than a conjecture at the 
moment, but is nonetheless suggestive. 

Now we turn to the precarious 0“^ theory for d = 3, and consider equations 
fl4.2.9p . (14.2.lip . (14.2.lOp and (I4.2.12p . For the theory is stable and 

asymptotically free and equation (I4.2.9P indicates that for a < 0, a is monotoni¬ 
cally increasing for fixed p as we increase the magnitude of the negative coupling 
(i.e., make the coupling more negative). Similarly for m? > equation (14.2.111) 
indicates that for a > 0, a is monotonically increasing for fixed p as we increase 
the magnitude of the positive coupling. For m? < the theory is unstable, 
and interestingly equation (14.2.lOp indicates that for such an unstable theory when 
a < 0, cr is monotonically decreasing for fixed p as we increase the magnitude of 
the negative coupling. Likewise for the unstable theory, equation (14.2.121) tell us 
that for « > 0 we have that a is monotonically increasing for fixed p as we increase 
the magnitude of the positive coupling. Thus, it appears that the direction of the 
monotonicity of a is tied to the stability of the theory. While we understand that 
the variational approximation is accurate for small values of |a|, we might trust in 
the variational approximation and conjecture that the sign of the derivative of an 
exact entanglement entropy with respect to coupling depends on the stability of the 
QFT in question. 
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Having examined the first derivative (|^)^ to analyze the monotonicity of the 

entanglement entropy with respect to coupling, we now turn to to analyze 

the convexity of the entanglement entropy with respect to coupling for a fixed value 
of the expectation value of the scalar held. Using equations fl4.2.2p and fl4.2.6p . we 
hnd that the second derivative of the normalized entanglement entropy per unit area 
a is 




a 


2(m^ — {m? — 


^ a {1 + a (1 + a 

Again using the ranges of the relevant quantities noted above, we have 

'd^a 


da^ 




a 






a 






a 


da^ 


> 0 


> 0 


> 0 


< 0 


1 + 1 dimensions 

2 + 1 dimensions 


(4.2.14) 



(4.2.15) 


(4.2.16) 

all a 

(4.2.17) 

all a 

(4.2.18) 


Similar to equation fl4.2.13p . we note that 


_ U / d^AS \ _ U /d^S\ 
da^J^ A \ da^ J ^ A \da^J 


1 + 1 and 2 + 1 dimensions 

(4.2.19) 


and so is a hnite and regularization independent quantity for d = 1 and 

d = 2. By equations (14.2.15p and (14.2.16p . for d = 1 and d = 2 we have that a is 
convex in the coupling for hxed +. Since for d = 1 and d = 2, a is monotonically 
decreasing, as the coupling increases the entanglement decreases less and less rapidly. 


Now let us examine precarious (j)‘^ theory where d = 3. For this theory, 
is not hnite and is not proportional to ■ In fact, the divergent term in 

is proportional to log(^)- Since for d = 3 we have 


(S). 

(ff). 


(9^(m^) 

da^ 


{w? — 1 ) 

a m^(l + a log(m^))^ 


(2a log(m^)^ + 2 log(m^) + l) — l) 

(4.2.20) 
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for small e we see that has the same sign as • Therefore, to analyze 

the convexity of S for precarious 0“^ theory it is sufficient to analyze the convexity 
of a. 

By equation fl4.2.17p . in precarious 0“^ theory for the stable region m? > 
we have that a is convex in the coupling for hxed <p. However, by equation (14.2.181) . 
in the unstable region where we have that a is concave in the coupling 

for hxed <p. We conclude that in the variational approximation, the convexity of the 
entanglement entropy is linked with the stability of the theory. One might conjecture 
that this connection between convexity and stability holds for exact calculations of 
the entanglement entropy. 
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Normalized Entropy per Unit Area a 


Entanglement Entropy in 1+1 Dimensions 



Figure 9: Normalized entanglement entropy per unit area in 1+1 dimensions as a 
function of the normalized coupling a, for various values of the normalized scalar 
held +. Lower lines correspond to higher values of the scalar held. 
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Normalized Entropy per Unit Area a 


Entanglement Entropy in 1+1 Dimensions 



Scalar Field 

Figure 10: Normalized entanglement entropy per unit area in 1+1 dimensions as 
a function of the normalized scalar field +, for various values of the normalized 
coupling a. Lower lines correspond to stronger couplings. 
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Normalized Entropy per Unit Area a 



Figure 11: Normalized entanglement entropy per unit area in 2+1 dimensions as a 
function of the normalized coupling a, for various values of the normalized scalar 
held +. Lower lines correspond to higher values of the scalar held. 
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Normalized Entropy per Unit Area a 


Entanglement Entropy in 2+1 Dimensions 



Figure 12: Normalized entanglement entropy per unit area in 2+1 dimensions as 
a function of the normalized scalar field +, for various values of the normalized 
coupling a. Lower lines correspond to stronger couplings. 
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Normalized Entropy per Unit Area a 



Figure 13: Normalized entanglement entropy per unit area in 3+1 dimensions, as 
a function of the normalized coupling a < 0, for various values of the normalized 
scalar held (p. The upper branches (in red) correspond to unstable solutions of 
the gap equation which do not minimize the Gaussian Effective Potential; the lower 
branches (in black) correspond to stable or metastable solutions of the gap equation. 
Curves from left to right correspond to increasing values of the scalar held. 
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Normalized Entropy per Unit Area a 



Figure 14: Normalized entanglement entropy per unit area in 3+1 dimensions, as a 
function of the normalized dual coupling a > 0, for various values of the normalized 
scalar held <p. Again, the upper branches (in red) correspond to unstable solutions 
of the gap equation which do not minimize the Gaussian Effective Potential; the 
lower branches (in black) correspond to stable or metastable solutions of the gap 
equation. Curves from left to right correspond to decreasing values of the scalar 
held. 
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Normalized Entropy per Unit Area a 



Figure 15: Normalized entanglement entropy per unit area in 3+1 dimensions for 
precarious theory, as a function of the normalized scalar held <p, for various values 
of the normalized coupling a. Couplings in this plot take only negative values; the 
coupling strength decreases in magnitude (becomes less negative and approaches 
zero) from left to right, or from inner to outer curves. The upper branches (in red) 
correspond to unstable solutions of the gap equation which do not minimize the 
Gaussian Effective Potential; the lower branches (in black) correspond to stable or 
metastable solutions of the gap equation. 
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5 Renormalization Group Flow of Entanglement 
Entropy 


In this section we compute the renormalization group equation for the normalized 
entanglement entropy per unit area a = —and show that it is simply 
related to the usual beta function for the renormalized coupling. For 1+1 and 2+1 
dimensional theory, as discussed earlier, there is merely a finite reparameteriza¬ 
tion of the coupling, with no coupling renormalization depending on a renormaliza¬ 
tion mass scale /i, so the beta function vanishes. However, for precarious theory 
in 3+1 dimensions, there is a non-vanishing beta function. If we differentiate the 
coupling renormalization equation 02.3.221) with respect to the renormalization scale 
H for hxed bare coupling and hxed cutoff, then in terms of the normalized coupling 
a = Aij/967r^ the beta function is given by 


P(a) = n— = 2a^ (5.0.1) 

a/i 

Since for precarious theory we have a < 0, the beta function in equation 05.0.11) 
tells us that the theory is asymptotically free. In particular, 1/a r'U log+2. 

Using the simplihed forms of the gap equation and entanglement entropy for¬ 
mulas in equations 04.1.16P and 04.1.17p . we obtain the renormalization group equa¬ 
tion for the normalized entanglement entropy per unit area 


da (l - ?) ^ da (l - ?) 1 


(5.0.2) 


We will refer to this quantity as the entanglement entropy beta function, although 
any such quantity obviously depends on the geometry of the region that dehnes 
the entanglement entropy as well as the quantum held theory itself. In 1+1 and 
2+1 dimensional theory, since the beta function vanishes, the renormalization 
group equation for the normalized entanglement entropy per unit area is trivial. 
Simplifying equation fl5.0.2p in the case of precarious (p^ theory, we obtain 


da 

l^-r = ^— -2^ 


(5.0.3) 


Since m?/pf < 1, this entanglement entropy beta function for precarious 0^ theory 
is “asymptotically free”, meaning that < 0, but only in the stable region of 
precarious <p^ where m?/p^ > Equality p^ = 0 holds only in the vacuum where 
mf /p^ = 1. The sign is reversed for the unstable region, with singular behavior when 
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m^/ jj? = where precarious 0“^ theory goes from being stable to unstable. It will 
be interesting to see which properties of the entanglement entropy beta function 
for precarious (j)^ are shared by other theories with asymptotic freedom, such as the 
Gross-Neveu model and non-abelian gauge theories. 


6 Conclusions 

We have applied variational methods to approximate entanglement entropies for 
scalar 0^ theories. Within our approximation scheme we have found that in d = 1 
and d = 2 dimensions, the entanglement entropy is monotonically decreasing and 
convex with respect to the coupling. This result is unexpected since the authors 
anticipated that larger coupling would lead to more entanglement. Our results are 
accurate for small coupling, although the approximation scheme is well-dehned for 
all values of the coupling. Inspired by the variational results, we conjecture that 
the exact entanglement entropy in d = 1 and d = 2 is monotonically decreasing and 
convex with respect the coupling for the entire range of positive couplings. Indeed, 
it would be strange if the derivative of the entanglement entropy switched signs for 
some particular value of the coupling. If this were the case, then this value of the 
coupling would be of immense physical interest since it would correspond to a local 
minimum of entanglement entropy. 

For precarious 0“^ theory in the stable range, we have also demonstrated that 
the variational approximation to the entanglement entropy is monotonic in a for 
a < 0 and for the dual values a > 0, and is convex for all a. This theory is of 
particular interest since it is asymptotically free. As we tune the parameters of the 
theory so that it becomes unstable, the derivative of the variational approximation 
to the entanglement entropy changes sign exactly at the transition between quasi¬ 
stability and instability. Likewise, the second derivative changes sign. We have also 
shown the the “entanglement entropy beta function” for precarious (f)^ theory in 3-1-1 
dimensions is asymptotically free. 

It would be interesting to see whether or not monotonicity and convexity with 
respect to coupling are generic features of entanglement entropy for QFT’s, and 
to explore why certain theories appear to have entanglement entropies which are 
monotonically decreasing with respect to coupling. It would also be interesting to 
examine the entanglement entropy beta functions of asymptotically free theories, and 
see if the behavior of precarious 0^ theory is generic. In future work, we plan to apply 
our techniques to quantum held theories with fermions as well as bosons. In addition, 
we hope to explore the entanglement structure of gauge theories with variational 
techniques. This would be a valuable step towards analyzing entanglement entropies 
of quantum held theories appearing in the Standard Model. 
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